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Abstract 

The magnetic and electronic properties of metal phthalocyanines (MPc) and fluorinated metal 
phthalocyanines (FieMPc) are studied by means of spin density functional theory (SDFT). Several 
metals (M) such as Ca, all first d-row transition metals and Ag are investigated. By considering 
different open shell transition metals it is possible to tune the electronic properties of MPc, in par- 
ticular the electronic molecular gap and total magnetic moment. Besides assigning the structural 
and electronic properties of MPc and FieMPc, the vibrational modes analysis of the ScPc-ZnPc 
series have been studied and correlated to experimental measurements when available. 
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I. INTRODUCTION 



Magnetism or magnetic moments are rarely found in organic compounds, but metal 
phthalocyanine (MPc) molecules are notable exceptions [Ij. The MPc family represents 
metal-organic semiconductors that display high thermal and chemical stability [2j. These 
molecules have recently attracted large interest in applications such as in solar cells, opto- 
electronic devices, light-emitting diodes, thin-film transistors, and gas sensors [3H7]. The 
magnetic properties of CoPc are used in spintronic devices [8j. Furthermore, the MPc are 
simple magnetic molecules which can serve as prototypical systems to study the magnetic 
properties of the 3d transition metals embedded in an organic surrounding. 

Fluorinated metal-phthalocyanines (FieMPc) are also widely used in applications. In 
contrast to the MPc which usually are p-type like semiconductors, FieMPc are used as 
n-type semiconductors in photovoltaic devices and in light-emitting diodes. Organic field- 
eflFect transistors have been intensively studied as components for low price, large area, and 
flexible circuit applications. To build them, an n-type organic semiconductor is needed. 
Experimental investigations focus on FieCuPc [9, lOj and FieCoPc [ ITHTl] . They have 
been also considered as an alternative to inorganic field-effect transistors (FETs) in some 
specific circuit applications such as radio frequency identification cards (RFIDs) p!5l [T6]. 
electronic paper [17], sensors jTSl and switching devices for active matrix fiat panel 
displays (AMFPDs) [20]. 

Several experimental techniques as X-ray spectroscopy f2T1-[25j. neutron diffraction [211 - 
[25] . scanning tunneling microscopy (STM) [26l[27j, nuclear magnetic resonance [28], and 
photo-emission spectroscopy [29j have been used to characterize the structural and electronic 
properties of MPc. However, for further developments in the synthesis and production of 
novel materials and devices based on MPc, more accessible and common techniques such as 
IR and Raman spectroscopy are useful. For example, Raman scattering has already been 
employed to identify different metal phthalocyanines and to probe polymorphic changes 
(packing) [30j, as well as structural modifications [3l]. Therefore, a systematic study of 
vibrational modes for the series MPc will provide information on the resolution of spectral 
features obtained with ultraviolet photo-electron, x-ray photo-emission and photoabsortion 
spectroscopy measurements. Vibrational coupling effects have been found to have an impor- 
tant effect in the description of outer molecular orbitals of CuPc [24], as well as of valence 
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bands and core levels of PbPc films [32]. The study of vibrational properties can help to 
improve the assessment of other theoretical methods used to study metal-organic complexes 
[33] . and to correlate magnetic transitions to changes of the surrounding [34] . 

From a theoretical point of view, density functional theory (DFT) is a widely used method 
to determine the electronic structure of various organic molecules, including MPc [3711351136] . 
Bialek et al applied an all-electron full-potential linearized augmented plane-wave method 
to study the electronic structure and magnetic moments of NiPc, CuPc, FePc, and CoPc 
[371440] . Marom et al have tested a variety of exchange correlation functional for copper 
phthalocyanine comparing local functional such as LDA (local density approximation) 
[42] and GGA (generalized gradient approximation) [431 144] with the semiempirical hybrid 
functional B3LYP [45j, the non-empirical PBE p6], PEE hybrid (PBEh) [47j, and the 
screened HSE [48]. While all these functional describe the geometrical structure quite 
well, the local functional underestimate the binding of orbitals at the molecule center. 
Additionally, GGA+U calculations have been performed in Ref. [49] to take into account 
the strong correlations that could be relevant around the metal center. 

The present article is organized as follows. After a brief description of the computational 
method in Section II, we report in Section III our results for the magnetic and electronic 
structure of CaPc, of the sequence of 3d MPc (M = Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn), 
and of Ad AgPc. In addition, we compare the electronic properties of the fluorinated MPc 
(FieMPc) and MPc. In Section IV we present our results for the vibrational properties of 
MPc, and we analyze general trends. A brief summary is given in Section V. 

II. COMPUTATIONAL METHODS 

In the present work, we focus on studying the electronic structure near the Fermi energy, 
Ef^ and the magnetic structure by means of density functional theory. Since hybrid func- 
tional with a stronger amount of Fock exchange tend to overestimate magnetic moments 
and magnetic order, we apply the simplest exchange correlation functional, namely the LDA, 
as implemented in the SIESTA package ^Oj. The electronic structure and the vibrational 
modes are calculated in MPc and FigMPc by using this implementation. The calculations 
have been performed by considering a finite spin polarization to allow for the formation of 
local magnetic moments. The PBE method [46] is also employed in order to compare with 

3 



the density of states (DOS) evaluated with LDA. 

Wavefunctions in SIESTA are described by a local atomic orbit als basis set; we have 
used a double zeta basis set and an energy cutoff of 300 Ry. The pseudo potentials are 
norm-conserving in fully non-local form. Structural optimizations were performed using 
the conjugate gradient method until net force on every atom was smaller than 0.04 eV/A. 
Non-symmetrical constraints were employed. At the optimized geometries, force constants 
were evaluated and employed to calculate vibrational modes and frequencies through the 
dynamical matrix. As we are interested in the electronic structure and vibrational properties 
of single molecules, we restrict our calculations to the F point, and a large simulation box 
of 28 A X 27 A X 45 A was used, to avoid periodic images interactions. Further geome- 
try optimization and frequency calculations of selected systems were also carried out with 
B3LYP/6-31G** in the Gaussian 03 suit of packages, in order to identify, by comparison 
with previous results, the IR and Raman representations. 

III. MAGNETIC AND ELECTRONIC PROPERTIES OF METAL PHTHALO- 
CYANINE 

Metal phthalocyanine, MC32N8H16, is an aromatic molecule with the molecular structure 
shown in Fig. [TJ Hexadecafluorophthalocyaninatometal, FieMPc, is constructed by replac- 
ing the hydrogen atoms within MPc with fluor atoms. It is expected that the metal (M) is 
located at the center of the molecule, displaying a square planar symmetry D4/^ [51j. The 
structural, electronic and magnetic properties of MPc are controlled by the particular char- 
acteristics of the metal within the inner ring. Phenyl groups contribute to the stabilization 
of these molecular complexes and influence their charge distribution, as reported for the 
formation of ScPc2 ^21 15^. Nonetheless, they may be subject to oxidation as in the case of 
TiPc [53H55] and VPc [Ml [57]. 

To study the magnetic structure of the flrst row transition metals and also of Ca and of 
Ag within an isolated MPc, we flrst focus on the determination of their geometrical struc- 
ture, and on how our observations relate to our approximations and the different molecular 
components. 

In general the optimized structures show a planar geometry except for the ScPc and 
CaPc, where the metal is slightly out-of-plane while the rest of the molecule undergoes a 
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FIG. 1. (Color online) Structure under investigation: metal phthalocyanine (MPc) molecule. The 
atoms are labeled according to the chemical constituents, and the metal atom in the center is 
denoted by M. 

small bending towards the opposite direction. This effect becomes larger as the atomic 
radius of the metal increases, the Sc atom (with radius of 1.62 A) is displaced by 0.24 A 
from the molecular plane while Ca (with radius of 1.97 A) is shifted by 1.12 A. To confirm 
the observed trend we calculate the optimized geometry for YPc obtaining an out-of-plane 
displacement of 0.72 A. It can be seen from Fig. [2] that for ScPc the valence electron density 
(isosurface value = 0.15) appears distant from the metal center, beyond the Van der Waals 
radius of Sc, i.e., the metal departure from the molecular plane is caused by the repulsion 
between inner electron shells on the metal and on adjacent atoms. In the distorted geometry 
the metal is more exposed and the molecular liability increases. 

Geometrical structures of minimal energy may be infiuenced by a Jahn- Teller distortion, 
which tends to reduce the symmetry from to D2h- As the point group symmetry was 
not constrained through our geometry optimizations, strictly speaking we have obtained, for 
MPc planar geometries, the D2h symmetry with two M-N distances slightly diflFerent from 
each other, nevertheless they differ only by less than 1 %, this variation being in order of 
magnitude comparable to the precision of the energy minimization procedure. A detailed 
account of the Jahn- Teller distortion effect on the studied MPc series is beyond the scope of 
the present work. Interested readers are referred to the work of Marom and Kronik who 
emphasize the theoretical and experimental difficulties to describe precisely the electronic 
states of MnPc and FePc, the differences being due to their intricate potential energy surface. 
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FIG. 2. (Color online) Charge density for ScPc. Isosurface value = 0.15 

Table H] shows the bond lengths between the metal and the nitrogen Nl, as indicated in 
Fig. [T] When available, we compare to experimental data. The theoretical bond lengths 
differ from experimental lengths by less than 0.05 A, i.e., by 1 to 2 %. Decreasing the atomic 
covalent radius leads to a decrease in the bond length from Sc to Cr, as well as from Mn to 
Cu. In FieMPc, the M-Nl bond length is typically slightly larger than in MPc, as shown 
in Tab. [1} The influence of the larger radius of fluorine with respect to hydrogen, and its p 
electrons show up only when no metallic states are present near the Fermi level. 

The last row of Tab. [l] shows the magnetic moments of MPc (in the case of FigMPc 
molecules, the calculated magnetic moments are quite close to the related MPc). Along 
the sequence from M = Sc to Mn, an open rf-shell on the metal progressively becomes half- 
filled and the MPc magnetic moments closely follow the tendency on the metal center (Sc 
{s = 1/2, 3rfi), Ti {s = 1, 3rf2), Y {s = 3/2, 3rf^), Cr {s = 2, 3rf^), and Mn {s = 5/2, 3d^)). 
The magnetic moment found for MnPc, 5/x^, is larger than the previously reported value of 
3/iB [59], however, high magnetic moments within this series are also obtained using PBE. 
In the second half of the rf-shell, total spin values lie between and 2. In general, they 
are smaller than those predicted for the first half. The values found for this interval agree 
very well with earlier results [21 [60l [6T]. Therefore, the trend in the magnetic moment can 
be correlated with the electron pairing and the crystal field created by the molecule to the 
metallic atom, which reduces the covalent character. The atomic orbitals hybridize such that 
they originate a new mixing, altering the expected orbital ordering from the plain crystalline 
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field effect, when it goes only over the non- hybridized atomic orbit als. As a consequence 
there is an abrupt decrease of magnetic moment from Mn to Fe, and zero magnetic moments 
are found for NiPc before completing the rf-shell. CuPc (5 = 1/2, 3d^) and AgPc (5 = 1/2, 
4rf^) have a hole in their nearly full rf-shell. The fact that magnetic moments for FigMFc are 
similar to those for MFc suggests that molecular magnetic moment are strongly influenced 
by contributions to the total electronic density coming from the metal center, with less 
dependence on the molecular environment; this observation holds for Ag and even Ca, an 
5-block element. 

The observed tendency for M-Nl bond lengths results from two main factors. The first 
one is the metal atomic radius reduction along the period, which clearly dominates from 
Ca to V, and the second one is the increase of electronic repulsion for high spin values 
and degree of orbital occupancy; these effects overcome the first one for Cr. In the second 
half of the rf-row, the decrease in magnetic moment values is reflected in the M-Nl bond 
lengths variation from Fe to Cu. Deviation from such behavior in the case of ZnPc may be 
understood by electron accumulation on the metal when the rf-shell is full, and the resulting 
increase of electronic repulsion. 



TABLE I. M-Nl bond lengths (A) and magnetic moments (m (/x^)) for MFc and FieMPc 



M 




Ca 


Sc 


Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


Ag 




MPc 


2.322 


1.998 


1.974 


1.968 


1.971 


2.012 


1.912 
1.927'^ 


1.883 
1.912'' 


1.875 
1.83*^ 


1.937 
1.935<^ 


1.975 
1.980^ 


2.042 


Bond 
























1.954-^ 






FieMPc 


2.328 


1.992 


1.988 


1.974 


1.975 


2.021 


1.991 


1.888 


1.884 


1.947 


1.988 


2.044 


m 


MPc 


0.00 


0.99 


2.00 


3.00 


4.00 


4.80 


2.00 


1.00 


0.00 


1.00 


0.00 


0.95 




abed 
1 1 1 


, *^ anc 


. f from references [22], 


M- EH, 123] and [S2 


, respectively. 



To complement the analysis of the magnetic properties of the MFc, we discuss their spin 
resolved electronic structure in the following. Firstly, we consider the density of states (DOS) 
of CaFc and FieCaFc shown in Fig. [3| Here and in the following plots, the density of states 
of one of the two spin components is multiplied by —1. Additionally, a level broadening of 
molecular orbital energies with a spread of 0.1 eV is used to simulate thermal and anharmonic 
effects in the spectra. The energy is shifted with respect to the Fermi level Ep^ estimated 
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using the Fermi-Dirac function with an electronic temperature of 300 K. The energy gap of 
CaPc, estimated as the difference between its highest occupied molecular orbital (HOMO) 
and its lowest unoccupied molecular orbital (LUMO), is 1.3 eV. As it can be noted in Fig. [3| 
these orbitals have no metal contribution, and they are localized upon the Pc atoms. Both 
peaks are very similar to the fluorinated derivative, the main difference is the decrease of 
the energy gap by 0.2 eV. However, the shape of the HOMO-1 (first orbital with an energy 
below the HOMO energy), is quite different: it is very broad in FieCaPc, in contrast to 
CaPc where it is rather narrow, indicating several energy states close by. 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 -2 -1.5 -1 -0.5 0.5 1 1.5 2 

E-EpCeV) E-EpCeV) 



FIG. 3. (Color online) Density of states CaPc (left) and for FieCaPc (right). 



To visualize the orbital contributions to the electronic density, we calculate the energy 
resolved density distribution nLDOs(r) = J^^ dE N{E^ r), where N{E^ r) is the local density 
of states. The HOMO and LUMO charge density isosurfaces for CaPc and FieCaPc are 
shown in Fig. [4| confirming that there is no contribution from the metal atom to these 
orbitals. This is related to the marked propensity of the alkali metal to loose its valence 
electrons and to form ionic bonds, transferring the charge to the MPc backbone. Dominating 
contributions in the case of the HOMO are from carbon atoms whereas orbitals localized on 
nitrogen and carbon contribute jointly to the LUMO. Therefore this state is more delocalized 
over the entire inner ring. The symmetry obtained for the HOMO's density indicates that 
it is an aiu orbital while the LUMO's density concurs with an Cg orbital type, the relative 
height of their corresponding DOS peak (see Fig. [3]) indicates that Cg is a doubly degenerate 
level. 



8 



FIG. 4. (Color online) Charge density isosurface for CaPc: HOMO (left), LUMO (right). 

The charge density isosurfaces for FieCaPc have almost the same shape as CaPc. In 
addition, the fluor atoms have clearly visible contributions. 

Turning to the 3rf-transition metals, the DOS of ScPc and its fluorinated derivative are 
shown in Fig. [5] and Fig. [6) The states around —1.1 eV are very similar to the HOMO of 
CaPc. The Sc d majority states show small contributions just below and above the Fermi 
level due to the mixing of the split dxz and dyz with the eg ligand orbital which, as mentioned, 
in CaPc is a LUMO doubly degenerate state. Metal contributions to unoccupied states also 
appear above 1.2 eV, coming from dxy^ d^z^ dyz and rf^2. The PBE DOS of ScPc indicates 
that the d^z contribution to the HOMO increases while it decreases for unoccupied levels 
appearing at lower energies (down-shift). A wider band gap of 0.3 eV is obtained when 
PBE is used. Regarding the electronic structure of TiPc (Fig. [7|), it can be seen that the 
metal has more contribution to the HOMO through its dxz and dyz orbitals than in ScPc. 
PBE predicts a slight HOMO down-shift, a marked peak splitting of dxy and rf^2 metallic 
contributions and a band gap increase of 0.3 eV. These effects concur with the expected 
reduction of the self-interaction error for the semi-local method. 

The electronic structure of FieScPc (Fig.[6]) indicates that states near the Fermi energy are 
not influenced by fluorination. In contrast, the main difference between TiPc and FieTiPc 
(Fig. [s]) is the level splitting of states with metallic contribution, particularly near the Fermi 
energy. The change in the admixture of Ti and Pc states modifles the majority DOS proflle. 
Even though fluorination can change the electronic structure, the total magnetic moment is 
not altered. This is due to the charge redistribution in the molecule, as can be concluded 
from peak intensity changes of electronic structure states close to Ep. The charge deflciency 
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on the metal and its magnetization, estimated as the difference between electronic charge 
densities of opposite spins obtained with the Mulliken population analysis [65], are reported 
in Tab. [11} The table shows that, with the exception of Ti systems, the small amount of 
charge transferred from the metal to Pc is similar for the fluorinated and not fluorinated 
compounds. When the tendency of estimated metal magnetizations is compared to that 
of the total magnetizations, previously shown in Tab. [TJ it is possible to observe a close 
interdependence. 



Sc 

I -1.5 



-0.5 0.5 
E-Ep 




FIG. 5. (Color online) Density of states of ScPc with LDA (left) and FEE (right). 



As noted in the second row of Tab. [TT| the estimated metal magnetization difference 
between TiFc and FigTiFc is 0.23 fiB- Along with this increase there is a diminution of 
O.lle of the charge transferred from the metal to the rest of the molecule. This can be 
understood as a compensation mechanism that explains why both molecules possess equal 
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FIG. 6. (Color online) Density of states of FieScFc with LDA. 
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FIG. 7. (Color online) Density of states of TiPc with LDA (left) and PBE (right). 
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FIG. 8. (Color online) Density of states of FieTiPc with LDA. 

total magnetic moments, stressing the ability of these compounds to modify their spin 
electronic densities by means of metal-ligand charge redistribution. 




Considering the LDA majority DOS of VPc, see Fig. |9| we note that new metal orbitals 
contributing to the occupied states appear between the HOMO and the aiu molecular orbital. 
They correspond to rf^2, dxy and, to a lesser extend, to dx2-y2. The degenerate dxz and dyz 
metal orbitals in the HOMO of TiPc become split in VPc, and they compose its HOMO 
and LUMO. The band gap is reduced to 0.4 eV. In contrast PBE predicts this gap to be 
0.9 eV, displaying a down-shifted HOMO (up-shifted LUMO) with a larger (smaller) metal 
contribution. Peak hights of rf^2 and d^y are similar to the dyz peak at HOMO, hence they 
are relevant for electron filling along the d-row series. Their energy ordering and particularly 
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TABLE 11. Band gap, metal magnetization and metal charge densities of MPc and Fie MPc 



M 




Ca 


Sc 


Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


Ag 


Band gap (eV) 


MPc 
FieMPc 


1.5 
1.3 


0.2 
0.2 


0.9 
0.4 


0.4 
0.5 


1.3 
1.2 


0.14 
0.14 


0.4 
0.4 


1.2 
1.3 


1.5 
1.4 


1.1 
1.4 


1.5 
1.3 


0.4 
0.4 


Magnetization (hb) 


MPc 
FieMPc 


0.00 
0.00 


0.19 
0.20 


1.47 
1.70 


2.97 
2.98 


4.02 
4.03 


4.78 
4.83 


2.58 
2.57 


1.31 
1.30 


0.00 
0.00 


0.47 
0.47 


0.00 
0.00 


0.38 
0.38 


Charge density (a.u.) 


MPc 
FieMPc 


1.26 
1.27 


0.30 
0.28 


0.28 
0.17 


0.15 
0.16 


0.05 
0.07 


0.48 
0.48 


0.30 
0.31 


0.23 
0.24 


0.14 
0.15 


0.15 
0.16 


0.39 
0.39 


0.25 
0.26 



the position of dxy^ indicates that the crystalline field eflFect is not determinant. In the PEE 
description the rf^2 and dxy metal orbitals appear at lower energies of aiu- 

Following along the 3d sequence, the additional electron of Cr causes important changes 



in the electronic DOS, see Fig. [T0| Metal contributions to occupied levels appear below the 
HOMO. The dxz and dy^ orbitals which are split in VPc become degenerate as in TiPc. For 
the first time within this MPc series the closest orbitals to the Fermi energy show larger 
energy differences for the majority states than for the minority states. It can be seen that 
there are no metal contributions to HOMO and LUMO, they principally localize upon the 
inner ring. HOMO 1 and HOMO-2 (first and second orbital below the HOMO energy) show 
only a small down-shift with PEE but the band gap remains almost the same. 
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FIG. 9. (Color online) Density of states of VPc with LDA (left) and PBE (right). 



The densities of states of MnPc and FePc are shown in Fig. [TTJand Fig. [T2| respectively. 
Inside the —2.0 ... 2.0 eV window around the Fermi energy, the only Mn non-negligible 
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FIG. 10. (Color online) Density of states of CrPc with LDA (left) and PBE (right). 




FIG. 11. (Color online) Density of states of MnPc with LDA (left) and PBE (right). 



contribution corresponds to the dx2_y2 {pig) to the HOMO. Other metal contributions to 
occupied states, previously seen in VPc and CrPc, appear at energies about 0.5 eV below 
the HOMO-1 {^lu)' The LUMO displays only ligand contribution, then the closest orbitals 
to the Fermi energy have a metal-ligand character with a narrow band gap of 0.14 eV (within 
LDA), and 0.18 eV (within PBE). These values are close to the PBE result, 0.21 eV, of Stradi 
et al. [66j who used different double-zeta (6-31G(d,p)) and triple-zeta basis sets for localized 
Gaussian- type orbitals (GTO), as well as plane- wave basis sets. Their peak decomposition 
analysis of DOS indicates that the metal contributes to HOMO and LUMO through dy^ 
and dxz] this result holds for every method they tested although it significantly reduces for 
LUMO with PBEO and HSE06. Different functionals (PBE, B3LYP, PBEh and M06) have 
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FIG. 12. (Color online) Density of states of FePc with LDA (left) and PBE (right). 



predicted either of two configurations if symmetry is enforced, aig |(HOMO) Cg |(LUMO) 
and |(HOMO) eg |(LUMO) [58]. On the other hand, when symmetry restrictions 
are not used either of two configurations are obtained, eg |(HOMO) eg t(LUMO) and 
eg l(HOMO) eg |(LUMO), as reported in the same work. Qualitatively, our result agrees 
with those configurations with a metal (ligand) contribution to the HOMO (LUMO) state. 
Both theoretical works point out that the calculated eigenvalue spectra of isolated MnPc 
appears compressed (into a smaller energy window) by about 20 % with respect to the 
reported data of the UPS experiment for thin and thick films deposited on metal surfaces. 
A more direct comparison is not possible because data for the gas phase are not available, 
and our calculations are basically for isolated molecules. 

In FePc, the addition of one electron initiates the second half of the rf-shell, where now the 
atomic orbit als can have two electrons, bringing about more possibilities to accommodate 
electrons and increasing the electron correlation due to the electron pairing. The HOMO 
and LUMO can be identified as eg orbit als with dyz and dxz metal contributions, respectively, 
when using LDA, and dx2_y2 metal orbital and eg with PBE. Regarding the metal (ligand) 
character of HOMO (LUMO), our PBE and LDA result concurs with recent calculations 
[601 [671469] but deviates from another prediction [58j. Furthermore, several theoretical works 
devoted to the interpretation of experimental results diflFer in the assessment of metal or 
ligand character of HOMO and LUMO. The analysis of electron energy-loss spectral data 
for FePc films [70j strongly suggests that partially occupied states e^, with dyz and dyz metal 
contributions, lie in between aiu and eg ligand states, manifesting metal participation for 
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HOMO with an energy difference to LUMO {eg orbital on the hgand) of 0.45 eV, very close 
to our LDA band gap, 0.4 eV, and PBE, 0.5 eV, estimates. The occurrence of metal orbitals 
close to the Fermi level, particularly for HOMO, has also been concluded from the analysis 
and interpretation of photoemission and X-ray absorption spectroscopy measurements of 
FePc on Ag(lll) [71j. The metal-ligand character obtained with LDA for HOMO-LUMO is 
used to describe the charge transfer between FePc and Au(lll) [72] by means of the LDA+U 
method. Analogously the GGA+U method is employed to reproduce charge reorganization 
of FePc adsorbed on Ag(lOO) deduced from STM spectrospcopy [69j. Apart from the +U 
correction our results differ in a qualitative manner because we find rf^2 states close the 
Fermi level but not dxy^ instead we get d^^-y^. This may be due to a larger degree of 
metal-ligand hybridization for dxy. In contrast a ligand-ligand character is suggested to 
account for photoelectron spectroscopy experiments which help to determine [/gff within the 
GGA+U correction Nevertheless, recent total-energy density functional calculations 

together with density-matrix constraints indicate that isolated and columnar stacking FePc 
molecules (a-FePc) ground states have a substantial d metal orbital contribution to HOMO 

EHl. 



Figure [T3| displays the density of states for CoPc. The LDA band gap, 1.2 eV, is larger 
than the previously reported result, LO eV [8j, using LDA and GGA+vdW. The band gap 
obtained with GGA, L45 eV, is the same as the one predicted for a CoPc monolayer [40] by 
the full-potential linearized augmented plane wave (FLAPW) method, and very close to the 
GGA+U result of L4 eV [69j, but smaller than the calculated value of L96 eV [35j using 
the Vosko-Wilk-Nusair (VWN) potential [42j including exchange and correlation gradient 
corrections. In the majority component, LDA and PBE indicate that HOMO and LUMO are 



distributed over the Pc atoms. From the charge density isosurface (Fig. 14) of the occupied 
Co states from — 2 eV to we see that the contribution to the majority density principally 
derives from Co and N states while carbon atoms contribute to minority states. In regard to 
the LDA DOS, we find that the metal-localized occupied rf^2 orbital down-shifts in the GGA 
description, on the contrary the LUMO, d^y^ up-shifts and the 2e^ orbital becomes LUMO. 
The ligand {aig) ^ ligand (2e^) character of HOMO - LUMO, differs from the metal (rf^2) - 
ligand 00], ligand ^ metal ligand (le^) ^ ligand (2e^) [35] character of other works, but 
agrees with recent theoretical results [69] in the gas phase, carried out to account for STM 
experimental observations. 

15 




FIG. 13. (Color online) Density of states of CoPc with LDA (left) and PBE (right). 




FIG. 14. (Color online) Charge density isosurface for occupied states, red minority-spin and blue 
majority-spin electron density. 



From the NiPc DOS functions shown in Fig. [T5| it can be noted that there is no difference 
between majority and minority spin densities. The electron pairing that occurs in this close- 
shell system causes equal positions and the same orbital order for both spins. Metal orbitals 
forming occupied orbitals are found below the HOMO, while at the same time there are only 
very small Ni dxz,yz contributions to the LUMO orbitals giving a ligand - ligand character 
to HOMO - LUMO with the same band gap of 1.5 eV for both spins. This value is smaller 
than the one obtained with the Hartree-Fock (HF) method, 2.41 eV [38j, but larger than 
the LDA prediction, 0.7 eV, for the monolayer and crystal monoclinic phase in the same 
work. Our band gap is very close to the estimation with the VWN potential of 1.47 eV 
[35] . but that study also reported that the LUMO corresponds to a big orbital, like a later 
PBE evaluation [2]. Regarding the band gap value and LUMO+l's identity, our result 
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agrees with more recent GGA+U [69] and Spectroscopy-Oriented Configuration Interaction 
(SORCI) calculations [60j. 



The DOS functions of CuPc (Fig. 16) indicates that splitting between big t and big I 
increases when the PBE functional is used, making the aiu orbital the HOMO and leaving the 
big I as LUMO. This up- and down-shift of occupied and unoccupied orbitals, respectively, 
is consistent with the expected reduction of the self-interaction error for localized levels. 
The ordering agrees with recent GGA+U calculation including the 1.4 eV band gap [69j, 
and differs from SORCI results [60j, and previous VWN predictions [35j. 



h 


1 

4 


tot 
Ni 










\ 




~ f 


H 





-1.5 



-0.5 





E-Ep 



0.5 



1.5 




FIG. 15. (Color online) Density of states of NiPc with LDA (left) and PBE (right). 




FIG. 16. (Color online) Density of states of CuPc with LDA (left) and PBE (right). 



The last two structures we discuss are ZnPc and AgPc. For ZnPc, the Zn metal states 
overlap with C and N states at 1.2 eV for the majority and minority states, see Fig. [T7| The 
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C and N contributions to the states are observed above —1.2 eV, specifically at —0.75 eV for 
C orbitals and at 0.75 eV for C and N orbitals. The band gap is 1.5 eV, similar to the CaPc 
molecule, but smaller than a previously reported value [31]. For AgPc, the Ag states mix 
with the N and C states at —0.1 eV for the majority states, and at 0.7 eV for the minority 
states. The pure C states prevail at —1.3 eV, C and N states at 0.2 eV and above 1.4 eV; 



thus it is difficult to decide whether AgPc becomes metallic or semiconducting, see Fig. 18 
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FIG. 17. (Color online) Density of states of ZnPc with LDA (left) and PBE (right). 



The band gaps for FieMPc and for MPc are given in Tab. [TT} The gap of each FigMPc 
is less than or equal to the gap of the corresponding MPc, which can be correlated to the 
fact that fluor is more electronegative, thus the electronic density in the central region of 
the molecule is more strongly attracted to the surrounding shell. 
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FIG. 18. (Color online) DOS of AgPc. 
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IV. VIBRATIONAL PROPERTIES OF METAL PHTHALOCYANINES 



The comparison of the calculated frequencies with experimental data of ZnPc indicates a 
semi-quantitative agreement for the main characteristics of Raman and infrared (IR) spectra 



(Tab. III). The experimental Raman spectra presented in Fig. 7 from the work of Tackley et 
al. [30] displays three bands which are more intense than the rest; our results coincide with 
their positions within 10 cm~^ with the exception of the most intense band at 1505 cm~^. 



TABLE III. Frequencies (cm ^) of relevant vibrational modes of ZnPc, and comparison with CuPc 







ZnPc 




CuPc 




ExperimentaP 


Calculated" 


Calculated'' 


Calculated*^ 


Raman 












1505 


1517 


1566 


1609 




1447 


1438 


1448*^ 


1491 




1340 


1362 


1347*^ 


1395 




1144 


1004 




1351 






1134 




1169 






1278 




1035 




747 


735 


737 


763 


IR 












727 


703 


706 


730 




752 


741 




769 




780 


760 




786 



^From reference [30j, ^Our CA/DZP calculation 
^Our B3LYP/6-31G** calculation 
^Sorted only by numerical value. See text for explanation. 

Our calculated frequencies of 1448 cm~^ and 1347 cm~^ are close to the experimental 
bands at 1447 cm~^ and 1340 cm~^, however, their displacement eigenvectors have more 
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similarities with those found in previously reported vibrational modes at 1438 cm~^ and 
1278 cm~^, respectively (see supplementary information of [30]). The IR and Raman nature 
of each vibrational mode considered here is compared to and confirmed through the displace- 



ment vector analysis of the calculated results for CuPc, also shown in Tab. Ill, The match 
between each vibrational mode type (IR) between ZnPc and CuPc give us the confidence to 
extrapolate these results to the rest of the MPc studied here. 

TABLE IV. Frequencies (cm~^) of relevant vibrational modes for MPc 



Mode description 



Metal out-of-plane 
Asymm. breathing 
Symm. breathing 
Ip^ N-M-N bending 
Op^ asymm. M-N4 bending 
Op^ symm. M-N4 bending 



Sc Ti V Cr Mn Fe Co Ni Cu Zn 



85 122 108 144 54 110 168 184 148 110 

- 193 205 201 171 175 191 197 174 159 

- - 262 262 229 - 260 255 259 264 
194 203 209 210 202 208 243 243 220 218 

- 279 - 275 - 243 314 - 250 223 
342 351 381 345 - 361 372 394 331 226 



^ Ip = in-plane, ^ Op = out-of-plane 



TABLE V. Frequencies (cm ^) of relevant vibrational modes for MPc 



Mode description 



Sc 


Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


326 


341 


352 


357 


312 


309 


389 


345 


301 


241 


329 


342 


358 


363 


316 


363 


401 


351 


302 


245 


544 


545 


542 


546 


537 


511 


546 


542 


539 


533 


576 


574 


572 


576 


569 


575 


581 


578 


577 


576 



M-N2 stretching-bending 
M-N2 stretching-bending 
M-N4 asymm. stretching 
M-N4 symm. stretching 
M-N4 asymm. stretching 
M-N4 symm. stretching 



1325 1183 1205 1357 1303 1342 1359 1362 1351 1347 
1366 1337 1387 1401 1348 1375 1407 1407 1403 1409 
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(b) 




FIG. 19. (Color online) Selected vibrational modes of CuPc also observed in other MPc, and 
described as metal out-of-plane displacement (a), asymmetric (b) and symmetric (c) breathing, 
in-plane N-M-N bending (d), out-of-plane asymmetric (e) and symmetric (f) M-N4 bending. 



Tables [rV] and \V\ show calculated frequencies of the MPc complexes for relevant vibrational 
modes such as breathing, and for those which involve atomic displacement near the metal. 



Displacement vectors for such modes are schematically presented in Figs. [19] and 20, As 
expected, symmetric vibrational modes are found at higher frequencies than asymmetric 
modes, but their frequency shift, when the metal center is varied, depends on the particular 
mode. Symmetric and asymmetric breathing movements of the entire molecule show similar 
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(f) 



FIG. 20. (Color online) Selected vibrational modes of CuPc, also observed in other MPc. M-N2 
stretching-bending, (a) and (b); M-N4 asymmetric ((c), (e)) and symmetric ((d), (f)) stretching. 



magnitudes of frequency shifts, the main difference being that frequencies of asymmetric 
breathing modes increase when going from Zn to Co, while the opposite tendency is found 



for the symmetric modes (see Fig. 21). This behavior may be explained by the release of 



stress around the metal inside the ring when it opens simultaneously in all directions. 



In the case of the M-N4 stretching modes, it is apparent from Fig. 21 that there is 
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a consistent increase from Zn to Co or Ni for the asymmetric type. The high-frequency 
(> 1180 cm~^) modes of M-N4 stretching have larger shifts. Their frequencies decrease 
from Fe to Co, then they increases again for Cr. These variations are not present for low 
frequency (< 600 cm~^) M N4 stretching modes: this can be understood on the basis of 



their displacement vectors (Fig. 20), which indicate significant deformations of benzene rings 
and smaller displacements on atoms close to the metal. Vibrational modes which involve 
atomic displacements near the metal show larger shifts. For bendings the magnitude of the 
shifts can be tens of cm~^. There is a frequency increase from Zn to Co, then a decrease 
for Fe, and an increase again for Cr. This is found for out-of-plane M-N4 bending, and for 
M-N2 stretching-bending modes. 



• asymm. breathing 
■ symm. breathing 
T op. asymm. M-N4 bending 
4 op, symm. M-N4 bending 
X M-N2 stieching-beiKlitig 
■t- M-N2 strching-bendin^ 





M-N4 symm. streching 
M-N4 asymm. streching 
A—* M-N4 symm. streching 
•w—w M-N4 asymm. streching 



square root of atomic mass ( a.u. ) 



FIG. 21. (Color online) Frequencies of selected vibrational modes of MPc 



In Tab. |Vl] we compare our frequency values and the nature of vibrational modes of CuPc 
with those reported in |73j for the fiuorinated system FCuPc, particularly for those which 
were found to have a larger contribution to the energy transfer process. Frequency values 
for CuPc become similar to results for FCuPc when the vibrational modes involve a small 
number of atoms, and when the atomic displacements are mainly from atoms closer to the 
metal. The above trends are also obtained with PBE. 
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TABLE VL Comparison of frequencies (cm ^) of relevant vibrational modes of CuPc and FCuPc 



FCuPc-^ 


175 983 1445 1570 


CuPc'' 


201 1158 1401 1561 



^From reference [73j 
^Our CA/DZP calculation 

V. SUMMARY 

The magnetic and electronic properties of metal phthalocyanines (MPc) and fluorinated 
metal phthalocyanines (FieMPc) have been studied by spin density functional theory for the 
first row-transition metals. CaPc and AgPc were also considered for comparison. A planar 
geometry with the metal lying inside the inner ring was obtained for all MPc except for the 
cases M = Ca and Sc, which are displaced from the plane by 1.12 A and 0.24 A, respectively, 
in direct relation with their atomic sizes. The analysis of electronic states shows that the 
total magnetic moment is caused mainly by the electron density localized on the metal and 
its nearby atoms. Furthermore, when the hydrogen atoms are replaced by fluor atoms the 
magnetic moments do not change. The same behavior is found for the band gap values. The 
M-N bond lengths of FieMPc are slightly larger than the corresponding ones of MPc. 

On the other hand, the analysis of the density of states (DOS) indicates a substantial 
change in the number, position, width and height of the DOS peaks in FigMPc in com- 
parison to MPc. The fact that the total magnetic moment values do not change can be 
explained by the charge density redistribution between the metal and its nitrogen and car- 
bon neighbors. This rearrangement suggests a high flexibility of the electronic distribution 
which may account for the apparent independence of the MPc properties on the metal, also 
pointed out in previous works [74j. However, the DOS proflles obtained along the studied 
d-Tow series indicate diflFerent mixing between metal and ligand orbit als, and corresponding 
energy shifts. 

The crystal fleld approximation cannot account for the orbital mixing occurring in these 
compounds. The ligand orbit als aiu and Cg have the propensity to appear near the Fermi 
energy and they can become HOMO or LUMO when the electronic repulsion shifts localized 
states on the metal or when electron pairing on the metal compensates repulsion between 
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levels for both spins as in NiPc and ZnPc. Near the middle of the studied d series, repulsion 
among unpaired electron densities causes different DFT functional to disagree with respect 
to the electronic structure. However, we have found that in general, LDA results are valuable 
for the qualitative assessment of the metallic or ligand character of the HOMO and the 
LUMO. Comparison with PBE results is useful to determine which localized metal orbitals 
are more strongly affected by the self-interaction error. The observed tendencies may be 
extrapolated to understand and predict the device performance of particular MPc employed 
in organic solar cells in which the exciton diffusion is affected by the energetical ordering of 
the orbitals [60j. 

The band gap shows a steep drop from CaPc to ScPc, i.e., when the d shell is initated. 
The largest value found in the d'^-d^ interval is for CrPc, and another important decrease 
is predicted for MnPc. The semiconducting character is recovered for the d^-d^^ series with 
low spin configuration. The crucial effect of the central metal in the electronic properties of 
MPc helps to explain recent theoretical results about the energetics of their adsorption on an 
Au(lll) surface as well as the increasing impact of the 3d states along the ZnPc-MnPc 
series as determined by energy-loss and photoemission spectroscopy [75j . Ongoing research 
regarding electron transport on these MPc is being carried out in our group. 

Based on the vibrational mode analysis, we find that, in general, frequencies increase when 
going from Zn to Fe if they correspond to asymmetric modes. There is an increment from 
either Mn or Fe to Cr. The series V, Ti, and Sc shows a more irregular behavior. However, 
there is a consistent frequency decrease from V to Sc for vibrational modes which do not 
include breathing or metal out-of-plane movement. The vibrational frequencies change a few 
tens of cm~^ when the metal in the central position is varied. This is an order of magnitude 
larger than the shifts provoked by molecular adsorption [76] . 
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